Electronic Structure and Doping in BaFe2As2 and LiFeAs: Density Functional 

Calculations 
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We report density functional calculations of the electronic structure and Fermi surface of the 
BaFe2As2 and LiFeAs phases including doping via the virtual crystal approximation. The results 
show that contrary to a rigid band picture, the density of states at the Fermi energy is only weakly 
doping dependent and that the main effect of doping is a change in the relative sizes of the electron 
and hole Fermi surfaces as required by Luttinger's theory. This is a consequence of a change in As 
height with doping, in particular a shift of As towards Fe as holes are introduced in the Fe plane, 
as might be expected from simple ionic considerations. The main effect of doping is therefore a 
reduction in the degree of nesting of the Fermi surface. This provides a framework for understanding 
the approximate electron-hole symmetry in the phase diagrams of the Fe-As based superconductors. 

PACS numbers: 



I. INTRODUCTION 

The discovery of high temperature superconductivity 
in oxypnictide phases, prototype LaFeAs(0,F) by Kami- 
hara and co-workersi has sparked widespread interest 
in establishing the physical properties of these mate- 
rials and especially the mechanism for superconductiv- 
ity. This interest has led to the discovery of a num- 
ber of new phases including other oxy-arsenides that 
when electron doped have critical temperatures exceed- 
ing 55K^ and hole doped superconducting phases includ- 
ing (K,Ba)Fe2As2, with Tc=38K^ and Lii_^FeAs with 
Tc of 18K.^ Thus there are three families discovered so 
far: (1) the oxy-arsenides, spacegroup PA/nmms^ with 
rare earth oxide layers separating the FeAs layers com- 
mon to all of these superconductors, (2) the body cen- 
tered tetragonal, lA/mmm ThCr2Si2 structure materials 
(BaFe2As2), where the FeAs layers are stacked so that 
the As atoms face each other and the Ba atoms sit in 
the resulting 8-fold coordinated square prismatic sites 
between them, and (3) P4/nmms LiFeAs, which is like 
the oxy-arsenide structure but with the rare earth oxide 
layers removed and replaced by Li (see below), result- 
ing in a large decrease in the c-axis spacing of the FeAs 
layers. Stoichiometric (undoped) BaFe2As2 shows a spin 
density wave (SDW) type magnetic order and a struc- 
tural transition with similar ordering temperatures. This 
feature is common to most of these materials, includ- 
ing the oxy-arsenides^ However, there are also some 
strong differences between properties of the oxypnictides 
and BaFe2As2. Most notable is the fact that with the 
exception of one report, ^^ the oxy-arsenide superconduc- 
tors are generally electron doped, while (Ba,K)Fe2As2 
and Lii_2yFeAs are found to superconduct when hole 
doped. However, in both BaFe2As2 and LaFeAsO su- 
perconductivity apparently is associated with the sup- 
pression of the SDW. This is potentially of importance 
because there has been much discussion of the association 
between magnetic quantum critical points and high tem- 
perature superconductivity and in addition the specific 



spin fluctuations associated with the nesting that pre- 
sumably drives the SDW have been studied as a possible 
pairing mechanismJ^ 

Electronic structure studies ^ii^ii^ii^ show that the 
Fermi surfaces and band structures of these phases are 
qualitatively rather similar. In particular, all the com- 
pounds show small compensating electron and hole Fermi 
surfaces, but high densities of states. The hole Fermi 
surfaces occur around the zone center, and are generally 
derived from heavier (lower velocity) bands than the elec- 
tron surfaces, which are around the zone corner (M point 
in a primitive tetragonal zone). Within band structure 
theories, moment formation is governed by a Stoner pa- 
rameter, /, which takes values of 0.7 eV to 0.9 eV for 
ions near the middle of the 3d series (note that the effec- 
tive / can reduced by hybridization). While magnetism 
may occur with lower values of the density of states, it 
must occur within a band picture if the Stoner criterion, 
N{Ef)I > 1, is met (the non- spin-polarized electronic 
structure becomes unstable against ferromagnetism in 
this case). Here N{Ef) is the density of states (DOS) at 
the Fermi energy on a per atom per spin basis. The cal- 
culated values of N{Ef) for the undoped FeAs materials 
put them at the borderline of itinerant magnetism, and 
in fact they are near both ferromagnetism, as required 
by the Stoner theory, and checkerboard antiferromag- 
netism, where nearest neighbor Fe atoms are oppositely 
polarized. The ground states, both experimentally^i^i^ 
and theoreticall y-*^ -'^i-'^^i-'^^i-'^^i-'^^i-'^^ have a magnetic structure 
that corresponds to condensation of an M point (zone 
corner) SDW. 

The electronic structures all show a strong increase 
in the electronic DOS as the Fermi energy is lowered 
into the heavier hole bands. This feature is quite ro- 
bust in this family and is also the case for the corre- 
sponding phosphide compounds.— ^2^ Importantly, this 
strong increase in the DOS poses a conundrum because 
it means that within a rigid band picture hole doped ma- 
terials would be strongly magnetic, and therefore would 
necessarily have very different properties, including su- 



perconductivity, from the electron doped materials. On 
the other hand, from an experimental point of view the 
generic phase diagram is much more symmetric, with an 
SDW and structural distortion near zero doping, with 
apparently non-magnetic, and superconducting states for 
both hole and electron doped materials. Here we investi- 
gate this issue using calculations for doped BaFe2As2 and 
LiFeAs within the virtual crystal approximation. This is 
an average potential approximation, which is beyond the 
rigid band method. We in fact find strong non-rigid band 
effects that leave only a weak doping level dependence to 
the N{Ef) and provide a plausible explanation for the 
generic phase diagram of this family of materials. 
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II. STRUCTURE AND METHOD 

The present calculations were performed within the 
local density approximation (LDA) using the gen- 
eral potential linearized augmented planewave (LAPW) 
method, ^^ similar to those reported previously for 
LaFeAsO.12 LAPW sphere radii of 1.8 ao, 2.2 ao, 2.1 
ao, and 2.1 ao were used for Li, Ba, Fe and As, respec- 
tively. The effect of doping was included by the virtual 
crystal approximation, by varying the nuclear charge of 
Ba or Li. Importantly, the internal parameters, including 
the As and Li positions, were relaxed independently for 
each doping level. 

We took the experimental tetragonal lattice param- 
eters of room temperature BaFe2As2 of Ref. Q, and 
relaxed the internal coordinate of As using LDA to- 
tal energy minimization. The calculated coordinate is 
^As=0-342. This is noticeably lower than the reported 
value of zas =0.3545. In fact the As heights in these 
two structures differ by 0.16 A, which as discussed below 
is outside the normal range of density functional errors 
and is significant for the electronic structures. We also 
did calculations with the generalized gradient approxi- 
mation but find a similar discrepancy. Similar discrepan- 
cies have also been noted in the LaFeAsO series ^i2^ For 
LiFeAs, we used the experimental lattice parameters^! 
and relaxed the free internal coordinates. The experi- 
mental structure (spacegroup PA/nmms, 129) shows two 
Li sites: 2c, which according to Ref. [2J is fully occupied, 
and lies above and below the centers of the Fe squares 
opposite the As, as shown in Fig. [H and 2b which has a 
reported occupancy of 0.1 in Lii.iFeAs and lies directly 
above the Fe. We performed total energy calculations for 
stoichiometric LiFeAs with Li in the 2b and 2c sites and 
found that the 2c site is favored by 0.52 eV per formula 
unit. Accordingly, we performed our doping studies plac- 
ing Li in the 2c site. We note that while this is the oppo- 
site choice to that made by Nekrasov and co-workers^^, 
the electronic structure near the Fermi energy is practi- 
cally unaffected by this choice. In fact, when hole doped 
by Li vacancies one may expect Coulomb repulsions be- 
tween the like charged Li ions to favor partial occupations 
of the 2b site, though always with a higher occupancy on 



FIG. 1: (Color online) Crystal structure of BaFe2As2 (left) 
and LiFeAs (right) with the relaxed As and Li coordinates. 
Note that in both structures the FeAs layers are uniformly 
spaced along the c-axis direction. 



TABLE L Structures used in calculations for undoped 



BaFe2As2 and LiFeAs. 





BaFe2As2 


LiFeAs 


Spacegroup 


M/mmm 


PA/nmms 


a(A) 


3.9625 


3.776 


c(A) 


13.0168 


6.349 


ZAs 


0.3545 


0.2116 


ZLi 




0.6728 


d(As-Fe)(A) 


2.315 


2.403 



the 2c site, similar to what is found in Naa^Co02. In 
any case, the insensitivity of the electronic structure to 
the Li position supports the validity of the virtual crys- 
tal approximation. The structures that we used for the 
undoped compounds are summarized in Table [H 



III. UNDOPED REFERENCE SYSTEMS 

The calculated LDA band structures and electronic 
DOS for BaFe2As2 are given in Figs. [2] and [3l respec- 
tively. These were calculated using the LDA relaxed As 
height as given in Table [H The Fermi surfaces for the 
calculated As height are depicted in Figs. [6] and [3 As 
noted previously^'— They consist of electron cylinders 
at the zone center and hole cylinders and other sections 
that depend strongly on the As height at the zone cen- 
ter. For BaFe2As2 the hole Fermi surface flares out at 
kz=^/2 (Z), giving it a more three dimensional charac- 
ter than that of the Li compound. This flaring out is 
sensitive to the As height, and in particular is reduced 
as the As is moved away from the Fe plane. In the Li 
compound there are three hole sections. These are two 
very two dimensional hole cylinders and an inner surface 
that forms a capped cylinder centered at F. Calculations 
were also done for the experimental As height (Fig. [3|). 
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FIG. 2: Calculated LDA band structure of BaFe2As2 for the 
LDA value of zas- The lower panel shows a blow-up near 
Ef^ which is set to eV. The short V-Z direction is from 
[0,0,0] to [0,0,1/2] while the long T-Z is in the basal plane 
and runs from [0,0,0] to [1,0,1/2] in units of the body centered 
tetragonal primitive reciprocal lattice vectors. 




FIG. 3: (color online) Calculated electron DOS of BaFe2As2 
for the experimental (top) and LDA (bottom) values of zas, 
plotted on a per formula unit basis. The projections are onto 
the LAPW spheres. Since As has extended p orbitals that 
extend well beyond the LAPW sphere, this yields values that 
are proportional to the As contribution, but underestimate it. 



These lead to a significantly higher N[Ef)-, which will 
make the material more magnetic, as has been discussed 
for the oxy- arsenide materials J^i^ The band structure 
and DOS for LiFeAs are given in Figs. [H and [Sj respec- 
tively, again using the calculated internal coordinates. 
The calculated values of N{Ef) are 3.06 eV~^ per for- 
mula unit (two Fe atoms) both spins for BaFe2As2 and 
3.58 eV~^ for LiFeAs on a per unit cell (two Fe atoms, 
two formula units), both spins basis. For the experimen- 
tal ZAs=0.3545 of BaFe2As2 we obtain N{Ef) = 4.59 
eV~^ per formula unit. These values are lower than the 
corresponding value for LaFeAsO (2.62 eV~^ per Fe both 
spins) ^ but are still high enough to place the materials 
near magnetism. 

For BaFe2As2 we find an antiferromagnetic state cor- 
responding to the SDW to be the most stable of the 
states considered. For the LDA value of zas we do not 
find any instability against either a ferromagnetic or a 
checkerboard (nearest neighbor) antiferromagnetic state 
in the LDA. The SDW state, which consists of lines of 
parallel spin Fe atoms in the FeAs planes has a moment 



defined by the integral inside the Fe LAPW sphere, ra- 
dius 2.1 ao, of 0.7 jiiB' As mentioned, N{Ef) is larger 
when the As height is raised to the reported experimen- 
tal value leading to a more magnetic state. In this case, 
within the LDA we find a very weak instability of the 
non-spin-polarized state to ferromagnetism (0.3 /i^/Fe, 
0.01 meV/Fe). a stronger instability for the checker- 
board antiferromagnetism (1.6 /i^/Fe, 41 meV/Fe) and 
the strongest instability to the SDW, which is the ground 
state (1.75 /i^/Fe, 92 meV/Fe). For the experimental 
value of zas the state with ferromagnetic c-axis stack- 
ing of the Fe in the SDW state has higher energy, and 
is therefore less stable than the state with antiferromag- 
netic stacking, by 3 meV/Fe, while these two states are 
degenerate to the precision of the calculation for the 
LDA As height. Thus the lowest energy magnetic struc- 
ture that we find is the SDW ordering in the Fe layers, 
stacked so that Fe atoms directly above each other in the 
c-direction are antiferromagnetically aligned. This is in 
accord with recent experimental results^ The high sen- 
sitivity of the Fe moment to the ordering and the fact that 
the SDW state is so much lower in energy than any other 
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FIG. 4: Calculated LDA band structure of LiFeAs for the 
LDA internal coordinates. 




FIG. 5: (color online) Calculated electron DOS of LiFeAs for 
the LDA internal coordinates, on a per formula unit basis. 
Note that the formula unit for this compound contains one 
Fe, while that for BaFe2As2 contains two. 





FIG. 6: (color online) LDA Fermi surface of BaFe2As2 for the 
LDA internal coordinates, shaded by band velocity. 



FIG. 7: (color online) LDA Fermi surface of LiFeAs for the 
LDA internal coordinates, shaded by band velocity. 



magnetic state strongly imply that the magnetism is of 
itinerant character with a spin density wave driven by 
the structure of the Fermi surface rather than local mo- 
ment physics associated with Heisenberg type exchange 
couplings. 

These undoped electronic structures are similar to 
those obtained previously and show the generic features 
of the calculated electronic structures of the other FeAs 
materials - small Fermi surfaces, with hole cylinders at 
the zone center and electron cylinders at the zone cor- 
ner, high N{Ef)^ and a strongly increasing DOS below 
Ep. The somewhat lower values of N{Ef) compared to 
LaFeAsO would place BaFe2As2 and LiFeAs further from 
ferromagnetism. This combined with the fact that fer- 
romagnetic spin fluctuations are highly pair breaking for 
singlet superconducting states may partially explain why 
hole doping is effective in these non-oxide materials, but 
the large increase in the DOS below Ef works against 
this explanation taken by itself. 



IV. DOPED MATERIALS: VIRTUAL CRYSTAL 
CALCULATIONS 

In this section we report virtual crystal calculations 
for the two compounds. This approach differs from rigid 
band in that it includes the self-consistent rearrangement 
of the charge density, which is done in an average poten- 
tial. This approximation is justified since the Ba and 
Li states do not contribute in any significant way to the 
band structure near Ef and because the scattering due to 
disorder on the Ba and Li sites is expected to be weak in 
view of the fact that wavefunctions of the states near Ef 
are primarily Fe derived and therefore separated from the 
Ba/Li atoms. This is supported by the fact that, as men- 
tioned, we find no noticeable difference in the electronic 
structure as calculated with the Li atoms of LiFeAs in 
sites 2b and 2c. Importantly, in the present context, we 
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FIG. 8: Calculated As height as a function of virtual crystal 
doping, S in carriers per Fe. Negative S corresponds to hole 
doping. 
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FIG. 9: Rigid band DOS as compared to the virtual crystal 
N(Ef) as a function of doping level for BaFe2As2. Note that 
the doping levels given in carriers per Fe with negative in- 
dicating hole doping, and the N{Ef) are given on a per Fe 
atom basis. 



are able to relax the internal coordinates and in particu- 
lar the As heights in the virtual crystal approximation. 

The calculated As height above the Fe plane as a func- 
tion of doping is shown in Fig. [51 As may be seen, the 
As drops towards the Fe plane as holes are introduced. 
This is as might be anticipated for a primarily ionic situ- 
ation, since a more highly positively charged Fe layer will 
be more attractive to As^~ and since the effective size of 
transition element ions in solids decreases as the valence 
increases ^i2U The essential point is that in these materi- 
als the carrier density is relatively low with an electronic 
structure near Ep formed from band edges. As a result 
the DOS and Fermi surface are sensitive to small changes, 
including the As position. 

Figs. [9l and [TOl show the calculated virtual crystal val- 
ues of N{Ef) for various doping levels as a function of 
doping level. As may be seen the behavior of the virtual 
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FIG. 10: Rigid band DOS as compared to the virtual crystal 
N(Ef) as a function of doping level for LiFeAs. Note that the 
doping levels given in carriers per Fe with negative indicating 
hole doping, and the N{Ef) are given on a per Fe atom basis. 



crystal is very different from what may be expected from 
a rigid band viewpoint. In particular, N{Ef) is found to 
be very nearly constant over the doping range studied, 
while as mentioned, the rigid band point of view would 
lead one to expect an increasing N{Ef) upon hole dop- 
ing. Our result differs from that of a recent calculation of 
Shein and Ivanovskii^ who studied a supercell of compo- 
sition Bao.5Ko.5Fe2As2, and found a significant increase 
in N{Ef) relative to the undoped compound. We ascribe 
this difference to the fact that they used the fixed exper- 
imental As position for the undoped compound and did 
not relax the internal coordinates. In any case, we find 
that contrary to the rigid band picture hole doping does 
not bring the system closer to Stoner ferromagnetism, 
which as mentioned would be detrimental to supercon- 
ductivity. Thus the main effects on the electronic struc- 
ture are more subtle. 

Figs. [11] and [12] show the evolution of the non-spin- 
polarized virtual crystal Fermi surface for the two com- 
pounds. As may be seen the basic structure of the Fermi 
surface with hole sections around the zone center and 
electron sections at the zone corner is maintained for all 
the doping levels shown. The difference between the vol- 
ume of the hole and electron sections is required to cor- 
respond to the doping level by the Luttinger theorem. 
Also it may be seen that the size of the electron sections 
changes at least as much as that of the hole sections (note 
that in Fig. [11] the top view emphasizes the flaring out of 
the hole Fermi surface at kz=0.5, while the main cylin- 
drical sections are smaller). This is contrary to what 
might be anticipated considering that the electron bands 
are more dispersive than the hole bands and therefore 
should change less and again shows the importance of 
non-rigid band effects. 

The condition for an SDW instability of the param- 
agnetic state is that the real part of the susceptibility, 
x(q) should diverge at the nesting vector. This function 
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FIG. 11: (color online) Virtual crystal Fermi surface of 
BaFe2As2 as a function of doping level b using the relaxed 
As position for each value of b. The shading is by velocity. 
The left panels show a view along the kz directions, while the 
right panels are along kx- 



is given in terms of an integral over the Fermi surface 
for the bare Lindhard susceptibility, Xo(q), with an en- 
hancement factor that emphasizes peaks, much like the 
Stoner enhancement of the Pauli susceptibility for the 
ferromagnetic case. The interband Lindhard susceptibil- 
ity for an electron cylinder at the zone corner and a same 
sized hole cylinder at the zone center is a function peaked 
at the nesting vector, while if the cylinders differ in ra- 
dius by (5g, e.g. due to doping or extra hole sections, the 
function will have a flat cylindrical plateau centered at 
the former nesting vector with diameter 2^q. Thus as the 
cylinders become mismatched in size the maximum value 
of Xo(q) (the value at the nesting vector) will decrease, 
while a region of high constant Xo(q) will develop within 
^q of the nesting vector. 

We calculated the Lindhard function, setting all ma- 
trix elements to be equal (i.e. constant matrix element 
approximation yielding arbitrary units). Although the 
Fermi surfaces are not perfect cylinders and there are 
multiple surfaces of different size, we do in fact find this 
behavior for LiFeAs. It is also the case for electron dop- 
ing in BaFe2As2 and to a lesser degree for hole doped 
BaFe2As2, where the peak in Xo(q) becomes broader and 



FIG. 12: (color online) Virtual crystal Fermi surface of LiFeAs 
as a function of doping level b using the relaxed atomic posi- 
tions for each value of 8. The shading is by velocity. The left 
panels show a view along the kz directions, while the right 
panels are along kx- The outermost hole cylinder is very two 
dimensional and therefore is difficult to see in the view along 
kz- Accordingly we indicate its position by an arrow. 



smaller but there is no clear plateau region. For both 
LiFeAs and BaFe2As2 we find a decrease in the peak 
height of the Lindhard function as we dope with holes. 
This explains why the SDW is destroyed by hole doping 
in these materials. 



V. DISCUSSION 

There has been much recent discussion of the relation- 
ship between the SDW instability and superconductivity. 
Within an Eliashberg formalism involving spin fluctua- 
tions what enters the calculation of the pairing function 
is a double integral over the Fermi surface (like that in 
the Lindhard function) with a kernel containing the pair- 
ing interaction, F(q) ~ x(q)- Thus in the case discussed 
above, the interband superconducting pairing between 
the electron and hole Fermi surfaces will depend on an 
integral over a range of q of size comparable to the Fermi 
surface size, while the SDW instability will depend on 
the peak value of x(q) (in a real case spin fluctuations 



away from the SDW ordering vector will compete with 
the SDW so even for a given maximum value of the bare 
Xo(q) a system with a narrower peak will order first). In 
any case, if spin fluctuations associated with the SDW 
are responsible for pairing, one expects superconductiv- 
ity to compete with the SDW and to occur on both sides 
of the region where the SDW is stable. What our results 
show is that the properties are much more symmetric be- 
tween hole and electron doping than would be expected 
from a rigid band picture. 
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